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BRAG MINIMIZATION FOR WINGS IN SUPERSONIC FLOW, 

WITH VARIOUS CONSTRAINTS 
By Max. A. Heaslet and Franklyn B. Fuller 

SUMMARY 



Che minimization of lnvlscid fluid drag le studied for thin 
aerodynamic shapes subject to imposed constraints on lift, pitching 
moment, base area, or volume. The problem is transformed to one of 
determining a two-dimensional potential flow satisfying either Laplace 1 s 
or Poisson* s equations with boundary values fixed by the imposed condi- 
tions. By mewyia of Kelvin's urir^rmnn energy theorem for harmonic fields, 
a method Is given for approximate drag minimization In the case of given 
lift. For supersonic-edged wings with straight trailing edges, perfect 
analogies are established between cases Involving lifting and nonlifting 
shapes. Particularly simple results are derived for a family of wings 
with curved leading edges with lift specified and center of pressure 
fixed at the 60-percent-chord position. General relations involving 
span load distribution a n d integrated loading along oblique cutting lines 
are derived. The minimum drag for other plan forms Is determined and, 
in the case of nonlifting wings, difficulties associated with unreal 
shapes are discussed. 


INTRODUCTION 


The calculation of supersonic drag of wings or bodies and the 
reduction of the minimization problem to one of determining a harmonic 
function of the lateral coordinates was reported by Nikolsky in refer- 
ence 1. Details of this method were not given, but a procedure leading 
to the same end that makes use of control surfaces which are everywhere 
inclined at the Mach angle to the streamwise direction was given by Ward 
in reference 2. Further work on the subject can be found in references 3 
through 6. The drag Is expressed by a surface integral over a surface 
that envelops upstream- facing Mach cones springing from the trailing 
edge of the wing or body. Drag minimization can then be re-expressed 
as a conventional Isoperimetric problem once the desired constraints, 
such as, lift, pitching moment, base area, or volume, are represented 
in terms of integrals over the same control surface. Various forms of 
these representations have been given previously; for convenience they 
are given again here in uniform notation. 
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The several possible variational problems are set up in the form 
of partial differential equations (Laplace* s _or Poisson’s) in terms of 
the lateral variables. The region in -which the solution is to hold is 
determined by the original wing (or body) , and the boun d a r y conditions 
to be applied are determined by the constraints of the particular prob- 
lem treated. Thus, although the problem has been reduced by one dimen- 
sion, solution is visually not easy because the shape of the regions in 
Which Laplace’s equation is to be solved is not generally suited to 
mapping techniques, wid the regions are, moreover, not simply connected. 
Numerical and approximate methods still apply, however, and particular 
f ami lies of solutions are readily found. 

An approximate method for the problem of minimum drag with given 
lift is developed, based on Kelvin’ s principle of minimum energy of 
potential flows. The results involve the apparent area (or mass) of the 
boundaries, so in reality require knowledge of potential flows about such 
shapes. In practice, the procedure offers a means for rationed, approxi- 
mation by enabling one to substitute known solutions for apparent area 
for the unknown ones of the problem at hand, and some examples are given. 

In any control- surface method of calculation of integrated forces, 
etc., a certain amount of detailed information is lost. Thus, in the 
present case, local surface shape is not an element of the problem, 
except insofar as it influences, say, the base area or the volume. 

Therefore it is not, strictly speaking, possible to give assurance that 
a real wing can achieve the minimum drag given by a solution to the 
variational problem. This can happen when the conditions of the problem 
allow a solution that results in a crossing of the stream surfaces that 
form the wing. An example of the occurrence of this is given; namely, 
the yawed elliptic wing of given volume, with zero base area. 

When the variational problems with the various constraints are 
listed, there is noted a similarity between the problems for given base 
area and for given lift, and those for given volume and for given pitch- 
ing moment. Conditions under which the constraints are interchangeable 
in each of these two cases can be found; the result being that, for cer- 
tain wings with supersonic edges, base area can be replaced by lift (or 
volume by pitching moment) without altering the drag of the wing. Such 
analog results have also been noted between optimum slender bodies and 
optimum lifting lines in subsonic flow. 

A particularly simple solution can be determined exactly for one 
family of plan forms. This family has straight trailing edges and 
supersonic-type leading edges that are formed by sections of a hyperbola 
asymptotic to Mach lines. The constraints specify lift and center of 
pressure (or pitching moment) . This family has two interesting limiting 
cases; at one extreme it becomes a sonic -edged triangle, and at the other 
a two-dimensional case. The triangular wing with sonic edges has been 
treated by Germain (ref. 4) and, by means of an approximation procedure, 
in reference 6, and comparisons with these results can be made. The span 
load distribution of these optimum wings is found directly from information 



NACA TN 1(227 


3 


derived in the analysis. It is also possible to determine, less directly, 
additional knowledge of aerodynamic loading integrated along oblique 
lines. So long as the obliquity of the lines is such that the component 
of stream velocity normal to them is supersonic, the variation of the 
integrated loading in the stream direction is obtainable. Hie so-called 
chord load distribution is a special case of these results. 


IMPORTANT SYMBOLS 


A 

Qd 

cl 

Cl ,C<2 


D 

f(y,z) 

h 

l 

L 

Mxq,0) 

q» 

m 

M 

n 


P 


base area of a wing or body 
drag coefficient 
lift coefficient 

curves bounding a region in an x «» const, plane 
(see sketch (b)) 

drag of a wing or body in a supersonic flow 

characteristic surface springing from a trailing edge 
(defined by x m f(y,z)) 

deviation of body from control surface at trailing edge 

streamwlse extent of wing or body 

lift 

integral of local loading along oblique line 

tangent of the angle of sweep or yaw 
Mach number in the free stream 

pitching moment, positive for a nose-up moment, taken about 
the line x = Xq, y = z ■= 0 

Inner normal to a plane curve 

inner normal to a surface 

pressure 


<3oo 


load coefficient (upper- surface pressure minus lower-surface 
pressure divided by free- stream, dynamic pressure) 
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<3oo 
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x,y,z 

*o 

Uoo 

V 

r t 


X,H,or,r 


Ho 

P 

9 

X 

V 2 


free -stream, dynamic pressure, i PocAio 2 
arc length 

region In an x = const, plane 
Cartesian coordinates 

x intercept of cutting line (see eq. (40)) 
speed of free stream, parallel to x axis 
volume of a wing or "body 

space curve (defined in sketch (a)) 

Lagrange multipliers (see eq. (10)) 

angle between oblique line and y axis (see eq. (40) ) 
density 

characteristic surfaces (see sketch (a)) 
velocity potential 

perturbation potential on surface x = f(y,z) (see eq. (2)) 

d 2 d 2 

two-dimensional Laplace operator, — - + — — • 

oy 2 oz 2 


ANALYSIS 

Expressions for the Drag and for Various Constraints 


Consistent with linearized supersonic theory, the boundary conditions 
over a wing or body are specified on a cylindrical control surface with 
elements parallel to the stream direction. The positive x axis is 
alined with the free- stream direction, and the perturbation velocity 
components u,v,w induced by the wing or body in the x,y,z directions 
are given by the gradients of the perturbation velocity potential cp(x,y,z) 
that is, by cp^, cpi^, cp . The differential equation governing the flow 
field is y 
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where 0 2 = The symbol 

denotes free -stream. Mach number and, In 
general, the subscript oo will "be used 
to denote free-stream conditions, for 
example, Ua, and are free-stream 

velocity and density. As shown in 
sketch (a), the vertices of the down- 
stream facing Mach cones trace out the 
leading edge of the control surface E 0 
and the vertices of upstream facing 
Mach cones trace out its trailing edge. 
The surface Z 0 is in this way- 
enclosed "by the envelopes E x and Eg 
of two families of cones. The surfaces 
E x and E 2 intersect along the space 
curve Ti. As shown in sketch (b) , the 
projection of Ti in a plane 
x = const, is the curve C x . Similarly 
the control surface E 0 projects onto 
the curve C 2 , and, the surface Es 
projects onto the two-dimensional 
region S. 

Let the surface Eg he given by 
x = f (y,z) and the value of the pertur- 
bation potential on E 2 be denoted by 
X(y,z), that is. 




X(y,z) = cp[f(y,z),y,z] (2) 


Sketch (b) 


Wave drag D (plus vortex drag when lift occurs) is then given by the 
expression (ref. 2) 


- ~ JJ' (Xy 2 +X z 2 ) cos(x,!T)dZ 

(3a) 

% 



(3h) 


S C 2 


where cos(x,N) is the direction cosine relative to the x axis of the 
inner normal IT to the surface ^ two-dimensional 

Laplacian operator (o /dy 2 ) + (a /dz 2 ) and 9X/5n is the gradient of 
the function X along the inner normal to a curve in the yz plane. 
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In the succeeding paragraphs, formulae for "base area, volume, lift, 
and pitching moment -will "be given. These results are applicable to any 
case In which the object of Interest deviates only slightly from a cylin- 
drical control surface . In reference 3, applications Involving both 
planar surfaces and quasi -cylindrical bodies of revolution were given, 
but in this report the examples will include planar cases exclusively. 

Consider a wing or body having thickness but without camber. If the 
thickness is zero at the leading edge, and idle projection on a plane 
x « const, of the deviation in area between the nose and tail of the body 
is A, it was shown In reference 3 that this area is given by 


_ 1 _ 

Uco 

JL_ 

Uoo 


JJ (Xyfy+X z f z ) CO s (x, N) az 

% 

s 2 

Jft v*x a, dz + i / 


C3.+C2 



(4a) 

g- 

(4b) 

t »as 

ckl 

(4c) 


When the aerodynamic shape is a wing, .the curve C 2 is a segment of 
the y axis . 

In reference 3 the volume enclosed by a nonlifting wing or body was 
expressed In the form 




where h(f,s) is the deviation from the control surface at the trailing 
edge. Equation (5) -thus gives the volume of aerodynamic shapes that are 
either open or closed at the base. In two general cases, however, the 
first integral on the right in equation (5) can be deleted: first, when 

the wing or body closes, that is, when h(f,s) » 0; second, when the 
trailing edge is normal to the stream direction, and the origin of coordi- 
nates is either on the trailing edge, or the base area is zero, that is. 


either f = 0 or 


J' h(f,s)ds 


0. In these cases the volume formula is 
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% 


V 




X(23*+V*f*)ay dz 



(6a) 


1 

2Uoo 


JJ 

8 


J ***<*. ± /*»!«. = 

S C3+C2 


(6b) 


The base area or volume cannot be affected by the fore and aft placement 
of the origin used in the analysis; proper positioning can, however, 
sometimes be vised to economize in algebraic manipulation. 

Consider next a wing or body of zero tfal ckness and carrying a load 
distribution by virtue of its camber and angl e of attack. The lift L 
and pit ching moment IM about the line x — xjn, z m 0 can be expressed 
in terms of integrals over the control surface S 2 (see sketch (a)) in a 
manner dir ectly* analogous to that used for the evaluation of base area 
and volume. The lift formula (ref. 2 ) follows from the application of 
linearized theory to the force relation 


f»- ff p^dS- ff pq ( q* ‘ N)dZ 

where 

? vectorial force on configuration 
df local velocity vector 
unit inner normal, to E 
p local pressure 

p local density 

The desired expression is 



( 7 a) 


(7b) 
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The general moment Integral is 


M - JJ p (4 "S )dZ + JJ p(^xf )< 3 S 

Sj_4£g 

■where r m is the vector distance between the moment center and an 
integration element on the control surface. When !M is positive for 
a nose-up pitching moment, linearized theory yields 

M - PoJJco JJ X V*(zf)ay dz + PcoLJoo J X p^ f) - XmCOs(z,n)^jds ( 8 ) 


When lift L is zero, or the pitching moment is calc ulat ed about the 
origin, one gets 


M 


8(zf) 

5n 


- PocAx> JJ x ^(zfjdy dz + PoJJoo J X 

S ° 2 

« PooUoo JJ Zf v®x dy dz + J 

S C1+C2 


ds 


(9a) 

(9b) 


Solutions of Variational Problems 


The problem of minimizing drag under the constraint of given base 
area, or volume . or lift, or pitching moment can be set up with the aid 
of equations (3), (4), (6), (7), and (9). It is only necessary to apply 
standard variational methods to any one of the expressions 


11 - D-M 

1 2 - BmV 

13 = D-cL 

14 ■ D-t]M 




> 


J 


(10) 
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where A,n,fiT,T are Lagrange multipliers. Che sign of \i is taken 
differently so that all quantities A,p.,cr,T can he identified as char- 
acteristic parameters in the so-called combined flow field. Thus, jj. is 
identified with pressure gradient in the combined field, a negative quan- 
tity. Any of these problems can be combined with another. For example, 
if it be required to find nrinlrniTm drag with given base area and lift, the 
quantity to be minimized would be D-AA-crL and so forth. The results 
found by applying the variational procedure to the expressions (10) will 
be given next. Each is a two-dimensional potential- flow problem in the 
lateral variables y,z. 

Given base area: 



< + ± f )“° 

in 

S 

bn V 

ic + f) - 0 

PccUoo ) 

on 

C 2 


X O 0 

on 

Ci 


(lla) 


D o - A 


(Ub) 


Given volume (zero base area) : 

*“(* + sfe; fa 


) P PcJJoo 


S ( x + f2 ) 


in S 


on C; 


(12a) 


X a 0 


on Ci 


D = - | V 


(12b) 


Given lift: 


- 0 in S 

(X + IWjz) =0 on C 2 

on 

X = 0 on Ci 


1 


(13a) 
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D - | L ■ (13b) 


Given pitching moment (zero lift, or moment center at origin): 


v®(x + 

UooTZf ) 0 0 

in 

S 



! (x+ 

UooTZf) « 0 

on 

C 2 

> 

(14a) 


X a 0 
Dai 

on 

Cl 


(14b) 


In each of the problems listed as equations (ll) through (l4) , the 
possibility of obtaining a solution depends primarily on the boundary 
curve Cx. Several cases where exact results could be obtained were dis- 
cussed In reference 3. However, recourse to approximate methods is 
usually indicated. One such approximation. In which the wing is distorted 
by Bntall amounts In order to achieve a solvable boundary curve Ci, is 
discussed in reference 6 . Germain, reference 4, has used development in 
series. Another method, for the particular case of finding mini nrum drag 
at given lift, is discussed in the next section. 


An Approximation of Minimum Drag Due to 
Lift by Energy Theorems 

From equations ( 3 ), ( 7 ), and (13b) , minimum drag for given lift takes 
the form 

D - 1 *■ as - - 1 I ds 

s ^ 

= Pw ^°° cr r X cos(z,n)ds a 2 L (15) 

2 ^Cz 2 


The first integral can be interpreted as the kinetic energy of a two- 
dimensional, incompressible fluid field (p^ = density) occupying the 
region S. Since the flow field is also lrrotational, one is therefore 
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prompted to use classical methods of approximation to estimate the drag. 
Lord Kelvin* s minimum theorem provides one such approach. It is stated 
hy Lamb (ref. T, p. 57) as follows: 

The irrotational motion of a liquid in a multiply- 
connected region has less kinetic energy than any 
other motion consistent with the same normal motion 
of the boundary and the same value of the total flux 
through each of the several Independent channels of 
the region. 

Consider, as in sketch (c), the minimum 
drag problem that arises when the wing trace 
lies well wi thin the curve Ci. The poten- 
tial field X(y,z) will be assumed laterally 
symmetric. The boundary condition along the 
wing trace calls for a constant value of 
dovnwash. On the exterior curve Cx the 
potential X is a constant so that the 
stream lines of the disturbance field are 
normal to C*. When Cx is far distant 
from the wing trace, the flow field corre- 
sponds to the case of a two-dimensional flat 
plate moving downward in a fluid field of 
infinite extent. When Ci is at a finite 
distance, an add i tional condition is imposed and, as a first approx! nation 
to the effect, one might assume that the flow in the vicinity of the plate 
Is similar to the flow near a plate moving downward, at a different and 
slightly increased velocity in an unrestrained field. The magnitude of 
the difference in velocities must be determined by the boundary condition 
an Cx* 



In order to use Kelvin’ s theorem It Is 
convenient to start with the flow about the 
plate moving downward at a velocity of mag- 
nitude, say, aUoo-k. Sketch (d) indicates 
the stream lines and the curve Cx Is shown 
dashed. The perturbation potential is 

- - (aUoo-k) <px(y, z ) 


where qpx(y,z) Is the perturbation potential 
for unit translational velocity. So long as 
the conditions in the flow are unchanged 
along the plate and at infinity, the kinetic 
energy of this irrotational field is a minimum, 
following three potential fields: 


w-'U^cr + k 



Sketch (d) 
Consider, now, the 


(a) The one in sketch (c) with X(y,z) satisfying equa- 
tions (13a). Set X = -trUorfp. 
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(b) A flow with uniform vertical velocity k within the 
boundary Ci. Its perturbation potential is ® 2 « kz. 

(c) The flow associated with a downward velocity -k* 
of the boundary C x . Its perturbation potential 
is <t>3 = -k’ cfs , where qfe corresponds to unit 
translational velocity. 

The sum of these three fields, see sketch (e), can be generated by 
imposing constraints along C x in the field of sketch (d) so long as k* 




is chosen consistent with the flow conditions at infinity. In order to 
meet this demand one uses the fact that, far distant, the induced field 
of a solid in pure translation behaves as a doublet (ref. 7, p. 165) . 

When the body is laterally symmetric, the axis of the doublet coincides 
with the axis of translation and the doublet or field strength is propor- 
tional to AfQ where A is the additional apparent area of the body 
and Q its geometric area. For ® x and ® 3 to behave alike, therefore, 
the relation 

(tfUoo-kjA^co = k' (Ag^Q) (16) 


must hold, where 

Affj fO0 apparent area of wing or plate in an infinite region 

Ag^oo apparent area of the boundary C x in rigid vertical motion 

Equation (l6) thus establishes one linear relation between k and k* , 
and it remains to establish one more. 

In the original field of sketch (d) no discontinuity in normal or 
tangential velocity occur b along C x whereas the combined field of 
sketch (e) has discontinuities in both quantities. These discontinuities 
are a measure of the added constraints imposed on the original, field. 
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Suppose, now, that k and k* are fixed so that the average value of the 
combined potential vanishes. Defining this average difference as 



one has 


kQ + 



cos(z,n)ds 


kQ + k* 




0 


or 


- * ,a b,o= 


Equation (17) furnishes a companion relation to 


equation (l6) . 


Kelvin's theorem gives 



(IT) 


where the Integration extends over all two-dimensional space and the 
right-hand member evaluates the contribution of X, ® s , and ® 3 . By means 
of Green* s theorem, one has 



Sn 


+ ® 


a ® 2 a ® 2 

2 dn dn 


+ 





ds 


Since ® 2 and X are harmonic functions, the final two terms In the right 
Integrand, are equal. The Inequality thus becomes 


(crUoo-k) 2 ^ a k' 2 % j00 + + Js^Q - 2kaU 00 A w ^ B 


where Ay -g is the apparent area of the wing -boundary region with the . 
potential/ X(y,z) . Rearrangement of terms yields 


Aw,oo ■%,»+« ft ' 

If one divides by (oUoo-k) £ Ay B and uses equations (l 6) and (17), 
the result is f * 
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1 

4*1,3 



_lJ\ . -jJ] s 17 i - 

Ab,o<^v ^oj L\Q 



The sign of the bracketed term on the left cannot be established In 
general but, if the inequality supplied by the minimum energy principle 
Is used merely as an approximation formula, the relation can be simplified 
to give 


_L_a 1 - 1 2 - 

*V f B Q -%,00+Q -*W,oo 

Since from equation (15) minimum drag is 


( 18 ) 


D “ f L " T U “ 2 ° faA W,B 


one has 


1w 3oo A W,B ^00 \8 -%,00+Q ^, 00 / 


( 19 ) 


R. T. Jones (ref. 8 ) has shown that the equality sign holds in this 
relation -when the curve Cx is an ellipse and the trace of the wing 
extends between the foci of Cx, a situation that always applies when the 
plan form of the wing is an ellipse. Further. teBts of the usefulness of 
equation (19) will be given in later sections. 


APPLICATIONS 


In this division some particular applications of the methods discussed 
above will be given. The first example, the elliptic wing, serves to 
illustrate a difficulty that can arise in the control- surface method for 
calculating drag. Next, same analog results are given, and, then, a 
family of wings with supersonic leading edges, having some interesting 
limiting forms, is treated. The concluding section shows how loading, 
integrated along oblique cutting lines, can be determined from general 
integral relations that apply to supersonic flow fields. 
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The Elliptic Wing of Given Volume 


A wing of elliptic plan form leads to a 
boundary curve Ci that is also an ellipse. 
Sketch (f) shows the wing, the boundary curves 
Cx and C 2 , and region S. Let Ci be given 
by 


B 2 C 2 


1 


The solution of the problem where minimum drag 
with given volume is sought is eased by the 
fact that by manipulation of equation (6b) , 
using -the first of equations (12a), the 
expression for volume can be put in the form 

7 ' ' £/(* + afe *“) * ** (20 > 

s 


h — 2c — H 



Sketch (f) 


Thus, the explicit determination of X from equations (12a) is avoided, 
and the complicated expression for f need never be used - only the fact 
that it vanishes on the outer boundary Ci. It is found that 


X + 


— - — f 2 

2poJJoo 


E P S B 2 C 2 Af zf A 

2 PooUoo B^C 2 \Jo C 2 ) 


( 21 ) 


satisfies equations (12a) . Substitution in equation (20) yields 

V *L_ B 2 C 2 S 

8qoo B^C 2 

where S is the area enclosed by Cij S = jtBC. How, with the aid of 
equation (12b) , the drag can be evaluated; 


D _ 4V 2 ( 1 l’N 

< 3oo V* 2 C 2 / 


(22a) 


Finally, in terms of the original wing parameters a,c (see sketch (f) ) 
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B 2 «= a 2 + c 2 /p 2 

a 2 - c 2 /p 2 

D 4V 2 3 

900 (jrac)c 2 (p2 +c 2/a2) 3 / 2 

(22b) 

which agrees with the result of Jones 
In reference 9* 

Next suppose that the e lli ptic 
wing Is yawed. Then, as shown. In 
reference 3, the bounding curve r*i 
lies in a plane x b m^y. Equa- 
tions (12a) can be solved by the same 
device as used just above in the 
unyawed case, and the result is 

_L + 1_ 

_D_ = 8Vf _1_ B g C g 

a ”“’® 0 e * 2 + s £( 3 + J) 

(23) 



Sketch (g) 


which reduces properly to equation (22a) 
when mi = 0. In this case, the expres- 


sions for mi, B, C in terms of the 
wing parameters a,c and the yaw angle if (see sketch (g)) are, if the 
wing eccentricity is e 2 = l-c 2 /a 2 - ■ • 



mi b (e g sin 2\|r)“ 1 ^-(T^p^sin^) 1/2 
B 2 /a g b (2P 2 )" 1 £t+ (l^-P^^sin 2 ^) 17 2 J 
C^/aa b B 2 /aa - (l-e g cos^jf) 

where 

T b p 2 (l-e 2 cos 2 t) + (l-e^in 2 ^) 


( 24 ) 
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In -this solution, the additional constraint of zero "base area is enforced, 
as it was in the unyawed case (eq.. (22b)) where the result checks that of 
Jones. The latter result applies to an elliptic wing that zero base 
area, and, in fact, is closed all along the trailing edge, being of bicon- 
vex section. Jones also deduced, from considerations based on Hayes* 
principle of equivalent source position and the resultant evaluation of 
wave drag, that the optimum yawed elliptic wing has a biconvex section. 

Let the wing in unyawed position have upper and lower surfaces given by 

z _ ± * (1 . . 2 ^ 

2 \ a2 c2 J 

where t is the maximum thickness. The volume of this wing is 


V 



If this wing is yawed through an angle if (see sketch (g) ) , its drag can 
be found by a single integration (ref. 10) : 


_D_ 

Qoo 



de 

I-*(0,t) 


with m = tan ifr and 


( 25 ) 


'(0,t) * — — (aa*-m a c 2 )-2Pm(a 2 -c 2 )cos &+-P 2 (m 2 a 2 *-c 2 ) cos 2 # 

1+m 2 L 


For the present purpose, it will be sufficient to consider the case where 
the m ino r axis 2c is small relative to the major axis 2a. 'Phpm 


i ■ -sc r — - iC stags; 

< 3 » s&aA J Q (i.mp cos 0 )* (i-m 2 3 2 ) 7 / 2 

(26) 

and this formula holds for mf 3 < 1 . For m |3 > 1, or when the wing is yawed 
outside the Mach cone, the drag is infini te. Equation (26) shows that 
the drag curve has a strong singularity as m3 -»■ 1. 


Now consider the limiting form of the drag given by equation (23) 
as e 2 -»l, or c/a-^ 0 . It is 


D _ kV s 2(l+m 2 ) g (2-m a 3 2 ) 

q °° 318,4 (2*-2m 2 p 2 -3m 4 p 4 ) (l-m 2 3 2 ) 1/2 


( 27 ) 



18 


NACA TN 4227 


Comparison of tMs result with that of equation (26) shows the latter to 
he the greater except at m ■ 0 where they are equal. Thus, two minimum 
drags have "been found, ostensibly for the identical constraints, and one 
is larger than the other. Thus it can he concluded that the larger drag 
results from a stricter constraint. This is cert ainly the case since 
pointwise t railing-edge closure is indeed more restrictive than zero base 
area, for the latter condition can be met by wings whose upper a nd lower 
surfaces cross to give an area distribution of hang in g sign at the 
trailing edge. 

The present method for minimum drag calculation has thus led to a 
difficulty. The closure of a trailing edge is not an easy condition to 
express in terms of the control- surface integrals, so it does not seem 
.feasible to force point-by-point closure by added constraints of the sort 
available. The element of ■uncertainty introduced by this situation, as 
to whether a determined minimum is achievable with a real wing, is cer- 
tainly undesirable. However, the question only arises in no nli ftin g cases, 
since a negative ordinate on a cambered surface causes no unreality. 


Analogous Variational Problems 


In the previous section entitled "Solutions of Variational Problems" 
were listed the differential equations and boundary conditions that are 
satisfied by the optimizing potential X(y,z) on the rear enveloping sur- 
face E s (see sketch (a)). If attention iB limited to wings with super- 
sonic leadi n g edges, and the lift-area and moment-volume problems are 
considered in pairs, a similarity that showB up can be exploited. By 
virtue of the supersonic leading edges, the variational problems can be 
analyzed in the upper half -plane z > 0, remembering, of course, that the 
flow and wing have vertical symmetry in the thickness case and that the 
flow field is antisymmetric in the lifting case. Consider, now, equa- 
tions (lla) and (13a). The curve Cg will be a portion of the axis z ■» 0 
and the two problems will yield the same X for z ^ 0 , and hence the same 

■minimum drag if 

vef - 0 and *? Z J Z=0 = 2(^7 (28) 


It is found that a straight supersonic trailing edge yields a rear 
enveloping' surface obeying equations ( 28 ). Let the t lulling edge have 
the equation 


x-ky = a ; k < p 

The Za surface is an inclined plane, «nd 
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f(y,z) = kyfc^-k 2 ) 1 /a z+a 

*]« - ± 0 2 -^) 1/2 

Equations (28) give 

A o ± 2 q 00 (p a -k a ) 1/2 <j 

Since drag is a quadratic function of X, and since, from equations (lit) 
and ( 13 b), equal drag occurs if 



it suffices merely to set 

L 3 ±2q 00 (p a -k a )“ 1/2 A (29) 


% 


Thus, if the ml n't mum drag of a supersonic-edged wing with a straight 
trailing edge parallel to x <=> ky of given base area is known, then it 
is also the minimum drag of the same wing considered as a lifting wing 
with lift determined by equation ( 29 ) . 

In this case, one needs to establish agreement between the relations 




n(2P a +f V®f) 
1* ( ffz ) z=o 


2qeoT(z ^f+2f z ) 

l 

2<3doT ( f + Zf z ) gtaO 


J 


(30) 


Equations ( 30 ) are satisfied by the same wing as in the previous case, 
except that k must be zero; that is, the trailing edge is normal to the 
free stream. In such a case 


and 




± 


2<3oqT 

0 


M = ± 7 

0 


(31) 
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Hals equation gives the pitching moment (at zero lift) of an optimum wing 
having the same optimum drag as the same wing carrying volume V with 
zero base area. The wings have supersonic edges and, a straight trailing 
edge normal to the free stream. 

Results of this sort have previously been noted between optimum drag 
configurations in subsonic lifting line theory and optimum slender bodies 
whose wave drag is given by the von Kerman formula. A summary of such 
results can be found in reference 11. 


A Family of Wings With Supersonic Edges 


Consider the family of wings whose plan forms are all portions of 
the hyperbola asymptotic to Mach lines through the point (-d,0,0). The 
equation of the leading edge is 



p2y2 0 2dx + x2 
and the trailing edge Is 
X a 2 

where the quantities d and l are shown 
In sketch (h). The root chord of the 
resulting wing is 2. If d-> 0, the wing 
becomes a triangle with sonic edges, and 
If d/2 » 1, the wing has very large span 
compared with its chord. 

The surface S 2 (sketch (a)) is 
composed of a pair of inclined planes 

x a f (y,z) = 2-0 |z| ( 32 a) 

and the boundary curve Ci (sketch (h) ) is 
made up of two parabolas 


PV* a (2d + 2)(2-20|z|) 


( 32 b) 


If mini-mum drag for fixed lift and center of pressure is sought for 
the wings of this family, the variation leads to the problem 


dn 


V 2 (X+U C0 TZf) a 0 
X+Uoo( C~XmT ) Z+U m T zf J a 0 

X o 0 


in 

S 

on 

c 2 

on 

Ci 


(33) 
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■where x^ is the coordinate of the center of pressure. A simple, exact 
solution of equations (33) follows directly if 


= d + f * 


in which case 


X a - 


Uoo£L 


16 [ l(2d + 2) ] 

In this event, the drag is given by 


375 -pr [p 2 y e +(ai+ 2)(2p|z|-Z)] 


( 3 *) 


< 3 co 4 VH 


-t 


d + | l 


[Z( 2 d+Z)] S/a 

1 

where L is tne given lift. Since the wing area is 


( 35 ) 


i (d+2) * 7 z( 2 d + 2) - d 2 coslf 1 2±1 
P L d _ 


the drag parameter is, written in temns of £ = d/2 


Pp 

3 Gl 2 


9 1+ 3 g 


40 (1+2C) 3/2 


[(3*5) ^25 - 



( 36 ) 


Sketch (i) shows the varia- 
tion of the drag with £. 
This latter parameter is, 
in geometrical terms, 

P 2 p q/Z, where p Q is the 
radius of curvature of the 
leading edge at the apex, 
and Z is root chord. 

As ^ varies from 0 to 00, 
the plan form, ranges from 
a sonic-edged triangle to 
a wing that has nearly a 
parabolic I eading edge. 

In the limit £-*- 0 , 
when a sonic-edged trian- 
gular wing results, the 
value of the drag param- 
eter given by equation (36) 
is 
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(37) 

This value is in agreement with the result of reference 4 for center of 
pressure at 60 percent of root chord. Also, the approximate result for 
given lift alone from reference 6 is quite close to that of equation (37), 
"being 0.223. This would indicate that the center of pressure is near 
the 60 -pare en t - chord position for given lift. In fact, from the results 
of reference 4, it is found that the center of pressure for given lift 
alone lies at 63 percent of the chord. The drag of the sonic-edged wing 
can also "be calculated "by the use of the approximation formula, equa- 
tion (19) . If the boundary Ca. is taken as an ellipse as a first crude 
approximation one gets 


On 

— % = 0.225 
PCl 2 


—z— = 0.211 

3Cl 2 

which differs from the result of reference 4 by about 5 percent. 

At the other extreme, £ — >», the wing plan form is very nearly a 
parabolic segment. The center of area of such a parabolic segment lies 
at 60 percent of the chord, which indicates that the loading is uniform 
over the wing. This is the correct result for minimum drag with given 
lift in two-dimensional flow, and the wing is Indeed becoming nearly a 
two-dimensional case as £ — >00. 


Triangular Wings With Subsonic Edges 



The approximation derived above as 
equation (19) can be used to predict drag 
of triangular wings with subsonic edges. 
The curve C3. for such wings is shown in 
sketch (j). It is made up of a parabolic 
center section faired into ellipses . The 
equation to the upper parabola is 

l a -2PZz = p 2 ^ 

and the ellipses are given by 

(y±m2/2) g + o x 

(Z/20) 2 (l-map 2 )(z/20) 2 " 


Sketch (J) 
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The curve Ci will be considered. as an ellipse in order to find its 
apparent area; 


Ab #00 = (If 2p) 2 Aj >00 - (1/23)** (Imp ) 2 

The other necessary quantities are 


Aw,*, = (z/2P) 2 aJ j00 = ( Z/23) 2 4rtm. 2 (3 2 


Q = (2/23)^* 


(l/2P) 2 2 ^l-n^g 2 cos -1 m0 + ^ m3(9-m 2 p 2 ) 


Then, from equation (19), 


gp 

PCt , 2 




The results of equation ( 38 ) are 
shown plotted in sketch (k) . Also 
shown is the drag curve for a flat- 
plate wing. When mp->0 the pre- 
dicted minlTmirn is exact to the 
order of linear theory; when mj3— 1 
the appro xima tions used introduce an 
unconservative error, as was found 
in the last section in connection 
with the sonic-edged wing. 



Sketch (k) 


Integral Relations for Loading 


The solution of the minimum drag problem pr o v ides a knowledge of the 
perturbation potential over the rear Mach surface ; the value of the 
perturbation potential at the trailing edge is therefore known an* the 
spanwise load distribution follows directly. In the family of plan forms 
with hyperbolic leading edges, for example, equation (34) yields 
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X(y,o) . - i ^ [p^-l(2d + J)] (39) 

16 Ez(2d+ z)] 3/2 

and the span loading remains parabolic for all- members of the family. 

It Is also possible to relate the chordvise distribution of load to 
information supplied by X(y,z) , and the following development will derive 
results of this nature. The first objective will be to determine the 
m ag ni tude of the load distribution integrated airmg an oblique line that 
lies in the plane of the wing. The Un a may be assumed to pass through 
the point (xo,0,0) and is to be inclined always at such an angle that 
the component of the free-stream velocity vector normal to it is super- 
sonic; it thus lies in the so-called supersonic zone of silence that 
extends laterally from the point (xo,0,0) . The method to be used is 
general in nature but, to avoid complication, attention will be directed 
principally toward the particular family of wings treated previously. 

In sketch ( 2) let OAB 
represent a supersonic -edged plan 
form (not necessarily a triangle) 
and let the trailing edge AB be 
x normal to the x axis. The cutting 
line CD is denoted 

x-y tan = Xq (40) 

Sketch ( 2) ■where is the angle between the 

line and the y axis and tan 0. 
Two regions of interest need to be distinguished. The first (Region I), 
which corresponds to the geometry shown in the sketch, corresponds to the 
range of n Q and Xq for which the cutting line does not intersect the 
trai ling edge. The second (Region II) arises when an intersection with 
the trailing edge does occur within the confines of the plan form. 

Region I .- Through the line CD pass two planes that are tangent 
to the double Mach cone with vertex at (xo,0,0). The forward and rear- 
ward inclined planes are, respectively, £3 and £4 and are represented 
by the linear equations 



x-y0 cob 0+z0 sin 9 » xq (4la) 

x-y0 cos 0-z0 sin 9 «■ Xq (4lb) 


where 


tan e 0 cos 9 

The notation used in equations (4l) is a common one in supersonic 
aerodynamics and was introduced first by W. D. Hayes in his treatment of 
equivalent source and doublet distributions (ref. 13) . 
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If Green's theorem for equation (l) is applied to an arbitrary 
region enclosed by a continuous surface E, the fundamental integral 
relation 



(te) 


results. The derivative 5<p/^v is the gradient of potential along the 
conormal v with direction cosines Vi, Vz f Vs that are related to the 
direction cosines n x , n a , ne of the inner normal, to the surface E by 
means of the relations 


-n^ 2 = AqVi , 


n 2 m ApV a , ns m AqVs 


(^ 3 ) 


Consider now the region enclosed by the surfaces E x , Eg, Eg^ and E 4 . 
Each can be generated as the envelope of a family of Mach cones. In every 
case, therefore, Ao = p, the conormals lie wlnug the surface, and the 
el e men t of area can be expressed as dE = dv ds where the element of 
arc ds is normal to dv and lies on the surface. Locally, the surfaces 
are cones and the orthogonal curves to the conormals lie always in an 
x si const, plane. The line CD may be written 


and the line EP is 


x =* Xp+py cob e 
2 ■ 0 


(«) 


x o ( z sin 0+Xo+Py cos 6) /(l+ Bin 0) 
z o (z-Xp-Py cos 0)/p(l+ sin 0) 


(^ 5 ) 


An integration with respect to v in equation (45) then gives, after 
noting that cp - 0 on E*, 


sin 0 Jq 


pD 

/ cp(xo+Py cob 0,y,O)dy 
vJC 


1+sin 0 /* Z sin GfXp+Py cos 0 Z-Xp-Py cos 0 _ 

sin 0 J E 1 h- B in 0 9 p (If sin 0) 


If this equation is differentiated with respect to Xp, the desired 
integral relation is obtained 
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pH nF 

/ u(xo+Py cob e,y,0)dy =. - / X z (y,z)dy 

Jq 23 J E 




In linearized theory, the local aerodynamic loading on a wing is given by 
4p/qoo ■ ^ u (zaOf)/°oo* If the Integrated loading along the line OD of 

equation ( 1 + 0 ) is denoted by L(xo, 0 ) /q^, equation ( 1 * 6 ) assumes the 
convenient form 


l(xo,e) 

<3oo 


— r P Xz<ay 

WcoJz 


(W 



c T 5 vT 7 

\ X 


E__ 

D i'N 

X N o 

x-/Sy 

* i 

^ 



ID "y 


Region II .- The cutting line now intersects the tr ailing edge. 
Equations (1+1) again represent the forward and rearward inclined planes 

through. CD and these planeB intersect 
1 7 the surface Eg along two lines. 

X s' ^ \ Sketch (m) indicates the geometry and shows } 

\ / nQ. \ above , the plane z » 0 and, below, the 

projection of the figure as viewed in the 
direction of the negative x axis. 

Equation (1+2) is again applied to the 
region enclosed by the surfaces Ex, Eg, 

Sq, and Z 4 , integration with respect to v 
is carried out and followed by differentia- 
tion with respect to Xg. The analogous 
Sketch (m) relations to equations (1*6) and (1+7) then 

appear in the form 


r D u(x a +- 3 y cos 6 ,y, 0 )dy - -if X z dy + A F\ z dy 
OC 2 3J E 23 Jq 


Qco 


- r\ay + r H x z ay 

3UocJe PUooJp 


W 

(^9) 



The above results can now be applied 
directly to the particular family of plan 
forms considered previously. For a fixed 
cutting an gle |i 0 , the range of Xq in 
Regions I and II is as shown in sketch (n) . 
The values Xo, Xx, and Xg are fixed, 
respectively, by the intersection of the 
wing’s axis of symmetry and three parallel 
cutting lines. The first of these lines 
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is tan g ent to the leading edge of the plan form, and the other two pass 
thr o ugh the intersection points of the 1 ending and trailing edges. One 

has Xo - -d(l-sin 0) , X x = 2-(2 a + 2Zd) x/2 cos 6, X 2 - 2*(2 a + 22d) 1/2 cos 0, 
and the following definitions 

Region I: Xq ■ xo »* Xx 

Region H: Xx ^ Xo ^ X 2 

When Xo ° Xx, Region I ceases to exist and this occurs when the angle n 0 
satisfies the inequality 


tan ji 0 > 


p(2 g +2d2) 1/a 

2+d 


Pram, equation (34) , 
1 




'00 


*z 


_ 3P g L 


8 Z 3 '' 2 (2+2d) x/2 


«nd the integration of equations (4-7) and (49) is Immediate once the 
intersection points D, E, F, and H are calculated. These points are : 


D: 


E: 


y - (2-Xo)/P cos 9 

(2 + 2d) cos 0 --|2(2+2d) (l+Bin 0) [xctf-d(l-sin 0) ]J- 


1/2 


p(l+sin 0) 


F; 


H: 


(2+2d)cos 0 + -j2(Z + 2d)(lrf-sin 0) [xQ+d( l-sin 0) ]J- 


p(l+sin 0) 


r *\ • 

(2 + 2d) cos 0-- 2( 2+ 2d) (l-sin 0) [xo+d(l+sin 0) ] » 
5 » — - J 

P (l-sin 0) 


1/2 


The Integrated loadings "become, finally. 
Region I ; 


L(xq,0) _ _3_ _L_ [xo+d(l-sin 0)] 1/2 
*3oo t/2. ^ 2 3/,2 (l+sin 0)^ 2 


( 50 a) 
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Region H: 


Mxq,9) 

<3oo 


3 L 1 

<3oo |S/2 


<• 

-(xp + 2d) ^ [xo+d(l-aln 0)] 1/2 [ 

[2(Z + 2d)] 1/2 cos 0 2(Lfsln 0) l/2 

to 


[Xo+a(lrfBln 8)] 1/g 
2(l-sin 0) 1/2 


(50b) 



In sketch (o), these loadings are drawn for the three members of the 
family corresponding to 5 ■ d/z t* 0, 1.5, 5, with p «= 1 and p D = 0°, 
15°, 30°, 45° in each case. The dimensionless quantities plotted are 


P(*b*,n) 


Lfxo^O)^ t *3oo^ 

<W L 


*0* 


Xo 

Z 


It is to be observed that when p ■ 0, 0 a jf/2 and Region II ceases 
to exist. In this case the integrated loading Is precisely the chord 
load distribution of the wing. One then gets the parabolic distribution 

L(xo,jt/2) ^ gL /VS 1 '" 

Qoo 2Zqco \ Z ) 


For all other cutting angles the load distribution falls to zero at the 
two end points. In the case of the triangular plan form (d=0), the 
streamwise gradient of L(xq,0) has a square-root singularity at Xq » O 
for all cutting angles. This, in turn, requires a square-root singularity 
in the slope of the wing at Xq = 0. 

Another integral relation between the loading and geometry of 
supersonic-edged wings can be found. If equation (42) is applied to the 
wing surface, the forward Mach surface, and the plane given by equa- 
tion (4la), one findB, for Region I, 
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J fiD r\D 

wCxo+Pyxcos 0,y x )ayi = 0 / u(xo+0y 1 cos 0,yi)dyi (51) 

C 

and a similar result can. be derived for Region H. Equation (51) can 
also be written 

L( , ^tfr " a wzf c w ^ x o + Pyi COB e »yi) a y3. (52) 

Thus the integrated loading across the wing is proportional to the average 
value of the wing slope along the same direction. This result is not 
confined to optimum wings but is of interest in connection with the 
problem of finding the surface shape to support the optimum load 
distribution . 
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